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Abstract 

We analyze the electromagnetic scattering of massive particles with 
and without spin and, using the techniques of effective field theory, we 
isolate the leading long distance effects beyond one photon exchange, 
both classical and quantum mechanical. Spin-independent and spin- 
dependent effects are isolated and shown to have a universal structure. 
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1 Introduction 



There has been a good deal of recent interest in higher order corrections to 
Coulomb scattering. In particular the one-photon exchange approximation, 
which has traditionally been used to analyze electron scattering has been 
shown to be inadequate when applied to the problem of isolating nucleon 
form factors via a Rosenbluth separation. Inclusion of two-photon exchange 
contributions have been found to be essential in resolving small discrepancies 
with the values of these same form factors as obtained from spin correlation 
measurements [Ij. A second arena where two-photon exchange effects are 
needed is in the analysis of transverse polarization asymmetry measurements 
in electron scattering. Such quantities vanish in the one-photon exchange 
approximation meaning that the sizable effects found experimentally must 
arise from two-photon effects [2]. 

Much has been written about such higher order photon processes and a 
number of groups have undertaken precision calculation of such effects [3]. It 
is not our purpose here to attempt such detailed calculations or to confront 
experimental data. Rather our goal is to use the methods of effective field 
theory in order to analyze the very longest range (smallest momentum trans- 
fer) contributions to the scattering process. These long range components 
are associated with pieces of the scattering amplitude which are nonana- 
lytic in the momentum transfer, and most of them are also singular in the 
limit of a vanishing momentum transfer (the exception being part of the the 
correction to the spin-spin coupling component where an extra factor of q 2 
arises). Some of these corrections are classical (ft- independent) and behave 



as l/v — q 2 while others are quantum mechanical (ft-dependent) and behave 

as log— q 2 , where q 2 is the invariant momentum transfer squared [3]. Below 
we shall examine both types of structures in the context of the electromag- 
netic scattering of two distinguishable massive particles of unit charge e with 
and without spin. In this case the lowest order interaction, which arises from 
one-photon exchange, is the simple Coulomb interaction 



where a = e 2 /47r is the fine structure constant. We find that two-photon 
exchange processes at threshold (v — > 0) yield corrections of the form 
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where Ac and Aq are the coefficients of the classical and quantum corrections 
respectively and are evaluated below. 

To see how such terms arise, in the next section we sketch our calculational 
techniques in the context of spin-independent scattering. This is a problem 
addressed nearly two decades ago by Feinberg and Sucher using dispersive 
methods [5]. Even earlier Iwasaki had studied the classical piece of this 
problem using standard noncovariant perturbation theory [BJ. Our quantum 
corrections are found to agree completely with those found by Feinberg and 
Sucher. However, our classical potential is at variance with that found both 
in [5] and [6j [7] . The origin of these differences is found in terms of differing 
contributions from the iterated piece of the lowest order potential, which must 
be subtracted from the scattering amplitude in order to produce a properly 
defined higher order potential [8]. Our work has also been motivated by more 
recent calculations in gravitational scattering where corrections to Newton's 
law are obtained [9]. 

In the following sections we extend these effective field theoretic methods 
to the problem of spin-dependent scattering and demonstrate that the results 
are universal, in that they can be written in terms of forms which are inde- 
pendent of spin. The calculation of spin-0 - spin-1/2 scattering reveals new 
structures of spin-orbit character whose universal form is also obtained when 
we consider spin-0 - spin-1 scattering. The next extension consists of spin- 
1/2 - spin-1/2 scattering where again we find the universal spin-independent 
and spin-orbit pieces as well as new (presumably universal) spin-spin cou- 
pling interactions. Our results are summarized in a short concluding chapter 
and the calculational details are found in the appendices. In Appendix IDI we 
give generalized results for arbitrary charges and g-factors of the scattered 
particles and argue for a multipole expansion like scheme which explains the 
universalities found. 

2 Spin-Independent Scattering 

We first set the generic framework for our study. We examine the electromag- 
netic scattering of two charged particles — particle a with mass m a , charge 
e and incoming four-momentum p\ and particle b with mass raj,, charge e 
and incoming four-momentum ^3. After undergoing scattering the final four- 
momenta of particle a is p 2 = Pi — 1 an d that of particle b is p 4 = p 3 + q — cf. 
Fig. [TJ Now we need to be more specific. 
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Figure 1: Basic kinematics of electromagnetic scattering. 



2.1 Spin-0 - Spin-0 Scattering 

We begin by examining the electromagnetic scattering of two spinless par- 
ticles. The electromagnetic interaction follows from making the minimal 
substitution in the Klein-Gordon Lagrangian density, yielding 

C = (iD^yiDy - m 2 <f)^<f) (3) 

where iD^ = id^ — eA^ is the covariant derivative, and leads to the one- and 
two-photon vertices 

r ^\p2,Pi) = -ie(p 2 +Pi),, 

Tj$(p2,Pi) = 2te\ u (4) 

Single photon exchange then leads to the familiar amplitude (in Feynman 
gauge) 
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with s — (pi + p 3 ) 2 the square of the center of mass energy. 

One way to define the nonrelativistic potential is as the Fourier transform 
of the nonrelativistic amplitude evaluated in the center of mass frame. We 
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Figure 2: One loop diagrams in electromagnetic scattering. 



will use a symmetric center of mass iram^j] with incoming momenta p\ = 
p + q/2 and p% = —p\ = —p — q/2 and with outgoing momenta f>2 = p — q/2 
and Pi = —p + q/2. Conservation of energy then requires p ■ q = so that 
Pi = p 2 + q 2 /4: for z = 1,2,3,4 and q 2 = —q 2 . In the nonrelativistic limit— 
q 2 ,p 2 <C m 2 — the amplitude reads 

« M a> ( g _ig f ! + _gL + . . .) + J^L f o + ( i + 2 m F +■■■)+■■■ 

(6) 

yielding the potential 

d 3 q 
(27) 3 



V«(r) = - / ^°M {1 \q)e^ 



r \ m a rrib J m a rrib \ 2m z m b * 



(7) 



The first component of Eq. ([7j) is recognized as the usual Coulomb potential 
(accompanied by small kinematic effects) while the second piece is a short 
range correction. 

Our purpose in this paper is to study the long distance corrections to 
this form which arise from the two-photon exchange diagrams shown in Fig. 
|2j This problem has been previously studied by Iwasaki using nonrelativistic 
perturbation theory [6] and by Feinberg and Sucher using dispersive methods 

1 These symmetric momentum labels of the center of mass frame are chosen so that 
the leading order coordinate space potential is real in the calculation of spin-0 - spin-1 
scattering presented below. 
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[5]. Our approach will be to use the methods of effective field theory, wherein 
we evaluate these diagrams by keeping only the leading nonanalytic structure 
in q 2 , since it is these pieces which lead to the long range corrections to the 
potential. This nonanalytic behavior is of two forms — 



i) terms in l/\/—q 2 which are ^-independent and therefore classical 

ii) terms in log— q 2 which are ^-dependent and therefore quantum me- 
chanical 

The former terms, when Fourier transformed lead to corrections to the non- 
relativistic potential of the form V c i ass i ca i{r) ~ 1/r 2 while the latter lead to 
V qU antum(r) ~ h/mr 3 corrections. For typical masses and separations the 
quantum mechanical forms are themselves numerically insignificant. How- 
ever, they are intriguing in that their origin appears to be associated with 
zitterbewegung. That is, classically we can define the potential by measuring 
the energy when two objects are separated by distance r. However, in the 
quantum mechanical case the distance between two objects is uncertain by 
an amount of order the Compton wavelength due to zero point motion — 
5r ~ H/m. This leads to the replacement 

l l l h 



which is the form found in our calculations. 

The calculational details are described in Appendix [A] Here we present 
only the results. Defining 

S = —j= and L = log —q 2 

we have, from diagrams (a)-(e) of Fig. [2] respectively 

°.Mg(?) = — (-2L) 
i^p Tn a rrib 

2 

M®(q) = ^—(2L + m a S) 
L^p Tn a rrib 
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where s = 
(m a + m b 



- (pi + P3) 2 is the square of the center of mass energy and s 
1 2 is its threshold value. Summing, we find the final result 
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We observe that in addition to the expected terms involving L and S there 
arises a piece of the second order amplitude which is imaginary. The origin 
of this imaginary piece is, of course, from the second Born approximation 
to the Coulomb potential, and reminds us that in order to define a proper 
correction to the first order Coulomb potential we must subtract off such 
terms. Before performing the necessary subtraction, we also point out that 
the imaginary part arising at this order of the amplitude is very singular 
in the nonrelativistic limit, even more than the leading order amplitude. It 
stems from an overall phase of the amplitude [19] and thus cannot contribute 
to any observable since observables are proportional to |7W | 2 . At our order we 
have calculated Ai = ^fA^ + fA (y2 ^ where an observable such as a differential 
cross section has a leading piece 0(a 2 ) and the corrections we have calculated 
contribute to order 0(a 3 ), but the imaginary part does not contribute to 
0(a 3 ). 

In order to subtract the second Born piece, we will work in the nonrela- 
tivistic limit and the center of mass frame — pi + p% = 0. We have then 



s - s = 1\ m? a + pi\hn% + p{ + 2p x - 2m a m b 



(10) 



and 



m a m b m r 

: r^j 



(11) 

s - s po 

where m r = m a m b / (m a + m b ) is the reduced mass and Po = \pi\, i = 1, 2, 3, 4. 
The transition amplitude then assumes the form 
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For the iteration we shall use the simple potential 

a 
r 

which reproduces the lowest order amplitude for spin-0 - spin-0 Coulomb 
scattering — Eq. ([7]) — in the nonrelativistic limit and which reads in momen- 
tum space 



0t/(1) 



Pf 
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where we identify Pi = p\ and pf = P2- The second Born term is then 
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where 



dH 
(2tt) 3 

G {0) (£) = 



v£\i-p f )GW(i)°v£\p l -t 



(2 



(15) 

— (16) 

is the free propagator. Note that in Eq. (floT) we take both the leading order 
potential as well as the total energies E(p ) and E(£) in the nonrelativistic 
limit. The remaining integration can be performed exactly, as discussed in 
Appendix [Cl by including a "photon mass" term A 2 as a regulator, yielding 



d 3 



(2*)' |ft + 



H 



L m r 

q 2 Po 



+ ie\£-pi\ 2 + \ 2 

(17) 



which reproduces the imaginary component of Adiot(<l)i as expected!] 

In order to produce a properly defined second order potential °V^ 2 \f), 
we must then subtract this second order Born term from the second order 
scattering amplitude Eq. (fl2|) . yielding the result 

d 3 q 



V^(f) 



-iq-r 



(27TY 



»M\ 2 JM- Q ^v [ c\q) 



2 We have omitted the IR singularity in the limit A — * since it does not contain 
nonanalytic dependence on q 2 . However, we note that it is present in the iteration as well 
as in the amplitude of the box diagram, and it is easily shown to cancel for the potential. 
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-tq-r 



d 3 q 

a 2 (m a + m b ) 
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m a m b 



-S(m a + m b ) + -L 
7a 2 h 



(18) 

lm a m b r' L <0Txm a m b r 6 
The result given in Eq. ffT8l) agrees with that previously given by Feinberg 
and Sucher for the quantum mechanical — ~ 1/r 3 — piece but disagrees for 
the classical — 1/r 2 — term. The classical contribution has also been calcu- 
lated by Iwasaki [6], who determined zero for this second order potential. 
The resolution of this issue was given by Sucher, who pointed out that the 
classical term depends upon the precise definition of the first order potential 
used in the iteration [8] and on whether one uses relativistic expressions in 
the iteration. Use of the simple lowest order form Eq. ( 1T31) within a nonrel- 
ativistic iteration yields our result for the iteration amplitude given in Eq. 
( TTTj) and is sufficient to remove the offending imaginary piece of the scatter- 
ing amplitude. However, if one uses relativistic expressions for the potential 
and the relativistic form of the energy in the iteration then alternate forms 
result with a different classical piece of the potential. Moreover, there are 
ambiguities in the leading order potential used for the iteration so that a 
unique definition of the second order potential does not exist [8]. 

For example, Feinberg and Sucher [5] calculate the iteration amplitude 



dH 


(Pf 
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E(p ) 
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) + ie 



(19) 



using fully relativistic expressions for the total energies in the denominator 
E(po) = E a (p ) + E b ( Po ) = s /ml +pl + ^ml+pl and E{£) = E a (£) + E b {£) 
and a potential including relativistic corrections 



Oy(l) 
V FS 



1 + 



P 




P 



(20) 



E a (p)E b (p) t V E a (p)E b (p) 

where we do not display the short distance part of the potential and with hats 
denoting operators whose ordering matters. Sucher calls this a "Feynman 
gauge inspired" potential whose operator ordering would be a natural choice 
when working in Feynman gauge [8] . One can evaluate the iteration integral 
by keeping the leading relativistic modifications of our previous results. Thus 
the Feinberg-Sucher potential — Eq. fl20l) — becomes 

Pi +Pf' 



Pf 



Oyd) 



FS 



1 + 



2m a m b 
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Similarly, the energy difference appearing in the propagator becomes 
E(p )-E(£) ~ 



m a + m b + -P- + lh) 



pi 



pi 



2m 
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We can now perform the Feinberg-Sucher iteration integral given in Eq. ([19 
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Ami 
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4m? 



(m a + m b ) 
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mi 



m a m b 
3-^ 
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m a m b 



m n 
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Am a m b 
m a m b 



m a + m b 



(23) 



which agrees precisely to this order with the exact relativistic evaluation 
given in Ref. [5]. 

Subtracting from °J^itot(o) we find then the second order potential 



d 3 q 
(2^)3 
d 3 g 



-iq-r 
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which is the form given by Feinberg and Sucher [5] . 

On the other hand Sucher [8] also discusses an alternative version of the 
one-photon exchange potential which is a natural choice when working in 
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Coulomb gauge. It includes relativistic expressions, and we refer the reader 
to [5] for the details of this "Coulomb gauge inspired" order a potential °Vgp. 
We note that on-shell the two potentials coincide 
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but when used off-shell in the iteration then the second order amplitude 
becomes 



°Ampg,($) 
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where fully relativistic expressions for the total energies in the denominator 
were used. Subtracting from °A4^ t (q) we find then the second order potential 
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This vanishing result for the classical piece of the 0(a 2 ) potential was also 
obtained by Iwasaki [6] and by Spruch [7]. Comparing Eqs. (|T8|) . (|24l) and 
( |27|) we note that the quantum mechanical contribution is invariant — only 
the classical term changes. As we have seen, even two fully relativistic itera- 
tion treatments do not yield the same answer for the classical piece, and the 
reason is that ambiguities arise when defining an 0(a) potential operator 
°V ( - 1 \ We have used an on-shell matrix element to find whereas the 

"Coulomb gauge inspired" leading order potential of [8j is defined from an 
off-shell matrix element. When we use these leading order potentials in an 
iteration where we integrate over all possible intermediate states, terms that 
vanish on-shell can contribute and yield differing results for the iteration 
amplitude and thus for the (D(a 2 ) potential. However, since the potential 
itself is not an observable this is not an issue. What is an observable is the 
transition amplitude. In each case we find the same result 



°Mtot(q) 



d re 



iq-r 



°V«(f)+°^(r1 



+ °Amp 4 (<f) 
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(28) 
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In this paper, since we are interested only in the threshold behavior of the 
transition amplitude, we shall utilize the simple Coulomb form for the po- 
tential and a nonrelativistic iteration, since this is sufficient to remove any 
pieces of the amplitude that would prevent us from writing down a well 
defined second order potential. 

Before moving on to details, it is useful also to point out the parallels 
between our calculational methods and those of the effective field theory 
NRQED [11] which has been set up to analyze non-relativistic bound states. 
The latter involves a systematic expansion in powers of the relative velocity v 
of the two particles. In this picture the one loop corrections to the amplitude 
involve terms of order l/t> 3 and higher. For example, the calculation of 
Manohar and Stewart^ evaluates corrections for the scattering of a particle- 
antiparticle pair up to order v° [12]. In coordinate space such terms include 
a combination of both short and long distance corrections. Our calculation 
involves a different sort of expansion looking only at the longest range terms 
in coordinate space, and it is not optimized for bound states since we do 
not use a power counting based on the virial theorem. Correspondingly, in 
momentum space we look for the nonanalytic components of the scattering 
amplitude. 

In their work in NRQED, Manohar and Stewart have performed a one 
loop matching calculation to O(v ) [12]. They calculate the full QED ampli- 
tude as well as the amplitude in their formulation of NRQED — vNRQED— 
which describes interactions of nonrelativistic fermions, ultra-soft photons 
and soft photons. For our discussion, the essential difference between vN- 
RQED and full QED is that potential photons have been integrated out and 
their effects are described by effective four-fermion operators in vNRQED. 
The coefficients in front of these four-fermion operators are the potential of 
Manohar and Stewart. 

In this paper we do not match onto a well-defined theory such as NRQED. 
Instead, our matching corresponds solely to the subtraction of the second 
Born iteration amplitude from the QED amplitude before Fourier transform- 
ing to a second order potential in coordinate space. Thus we regard our 

3 Manohar and Stewart (and collaborators) have performed a number of NRQED and 
NRQCD calculations of great phenomenological importance [13], such as the Lamb shift 
and hyperfine splitting or the ti production cross section near threshold. Much of their 
work includes impressive two loop results and involved applications of renormalization 
group methods. When we refer to Manohar and Stewart in our comparison here, we refer 
specifically to Ref. [12] which we found most suitable when comparing with our work. 
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potential as a nice way to display our resulting scattering amplitudes in co- 
ordinate space, but we emphasize that our main results are the long distance 
components of the scattering amplitude. Since we do not match onto a the- 
ory containing photons whereas the Manohar-Stewart potential arises in the 
matching to vNRQED which has soft and ultrasoft photons as degrees of 
freedom, we expect that the potentials differ. It is seen that our quantum 
corrections proportional to 1/r 3 stemming from log q 2 pieces of the amplitude 
are absent from Manohar and Stewart's potential because the exchange of 
soft photons in vNRQED yield the complete QED contribution of log q 2 terms 
so that their matching does not yield any quantum pieces of the potential 
Despite these differences, Manohar and Stewart's calculation [12J of spin- 
1/2 - spin- 1/2 scattering gives us the opportunity to compare our results for 
the scattering amplitudes and iterations with theirs. Note that one further 
difference is that our calculation deals with non-identical particle scattering, 
while that of Manohar and Stewart involves quark-quark pairs or quark- 
antiquark pairs of equal mass. Thus we do not have the exchange piece of the 
amplitude or the annihilation channel contributions given in [T2]. We verify 
that the classical spin-independent iteration amplitude in vNRQED found by 
these authors using on-shell matching — given by the sum of iterated terms 
involving V c x V c (including relativistic corrections to the free propagator) 
and V c x V r in their notation — yields for the classical piece 



and agrees with the iterated Feinberg-Sucher amplitude Eq. (1231) when we 
set m a = rrib — m - Note that Manohar and Stewart also emphasize the 
nonuniqueness of the classical potential, which they associate in part with 
the existence of off-shell terms of the lowest order potentials such as 



4 The alternative formulation of potential NRQED (pNRQED) of Pineda and Soto 
[T3] , [T5] differs from vNRQED in that it only contains ultra-soft photons as degrees of 
freedom whereas soft photons are integrated out. (Further differences between vNRQED 
and pNRQED are not important for our discussion here.) Therefore, the potential in 
pNRQED does exhibit the same quantum pieces of the potential as ours — cf. Eq. (2.17) 
in Ref. [15] for example. However, we note that Pineda and Soto use off-shell matching 
in Coulomb gauge yielding a vanishing classical piece of the potential. 



(29) 




4m 2 q 4 
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which vanish on-shell. Since we perform the iteration using on-shell potentials 
these terms do not contribute, but they can contribute in some forms such 
as those used by Iwasaki [6] and in the Coulomb inspired form of Sucher [7] . 

However, despite the agreement of many forms, the corrections which 
we examine are often higher order in the relative velocity than included in 
the Manohar- Stewart vNRQED exposition of [12]. Thus in the case of the 
spin-orbit term, the one loop corrections which we consider are order v in 
the Manohar-Stewart expansion and are therefore outside the quoted pieces 
of their potential. Likewise, the spin-spin correlation corrections which we 
retain are order v in vNRQED and again are not found in the Manohar- 
Stewart expressions of [12J. 

3 Spin-Dependent Scattering: Spin-Orbit In- 
teraction 

3.1 Spin-0 - Spin-1/2 

Having determined the form of the scattering amplitude and the resulting 
potential for the spinless scattering case we move on to the case of scattering 
of particles carrying spin. We begin with the scattering of a spinless particle 
a from a spin-1/2 particle b. From the Dirac Lagrangian density 

C = -00)0' JP - m)ij){x) (30) 

we determine the one- and two-photon vertices for a spin-1/2 particle to be 

T<$fa,P3) = (31) 

and the resulting transition amplitude at tree level is found to be 

where our spinors are normalized as u(p)u(p) = 1. Defining the spin vector 
as 

Sg = -ufafrrfufa) (33) 
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where 75 = — Z7°7 1 7 2 7 3 , we find the identit;yE| 



1 - 

Ami 



(P3 + P4) 
2lTLh 



mi 



(34) 



whereupon the nonanalytic part of the transition amplitude can be written 
in the form 



Ana 



m a m\ 



(35) 



Now we again give the nonrelativistic amplitude in the symmetric center of 
mass frame (p\ = —p 3 = p + q/2) where 



S - ™> (0, S b ) with S b = - x pxl 
u{pMp,)^ X b } X \-^S b .p*q 



(36) 
(37) 



and 



ta^&PiPzq 1 S s b ^ (m a + m b ) ( 1 + 



P~ 



2m a m b 



S b -pxq, (38) 



so that we find 

^M^(q) 



Ana 



6f 6 , + 2m b) 3 - - 



2m a m 2 b 



(39) 



whereby the lowest order potential becomes 

d 3 q 



(2tt)< 



*M (1) (q) fT** 



^ -X/Xi-^; —S b -pxV- 

r 3 2m a m b r 

- -X/Xi-— -5 2-L-S b 

r 3 r 6 2m a m b 



(40) 



where L = fxp is the angular momentum and r = f a — f& — the modification 
of the leading spin-independent potential has a spin-orbit character. 



5 Here e 0123 is taken to be +1. 
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When evaluating the one loop corrections we encounter an additional 
complication: The calculation contains two independent kinematic variables, 
the momentum transfer q 2 and s — So which is to leading order proportional 
to Pq (where p = \pi\, i — 1,2,3,4) in the center of mass frame. We find 
that our results differ if we perform an expansion first in s — Sq and then in q 2 
or vice versa. This ordering issue only occurs for the box diagram, diagram 
(d) of Fig. [21 where it stems from the reduction of vector and tensor box 
integrals. Their reduction in terms of scalar integrals involves the inversion of 
a matrix whose Gram determinant vanishes in the nonrelativistic threshold 
limit q 2 ,s — SQ — > 0. More precisely, the denominators or the vector and tensor 
box integrals (see Appendix |AJ involve a factor of (4p§ — q 2 ) when expanded 
in the nonrelativistic limit. Since q 2 = ^p\ sin 2 | with 9 the scattering angle, 
we notice that 4p^ > q 2 unless we consider backward scattering where 6 = n 
and where the scattering amplitude diverges. And since p^ originates from 
the relativistic structure s — sq, we therefore have to first expand our vector 
and tensor box integrals in q 2 and then in s — Sq. Evaluating the diagrams 
(a)-(e) of Fig. [2] and we find the results 
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Summing, we determine 
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Using the identity Eq. (134"1) and 
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Eq. ( 1421) becomes 
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Finally, working in the symmetric center of mass frame and taking the non- 
relativistic limit using Eqs. ( 1371) . ( l38l) and 
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we find 
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We note from Eq. ( 1451) that the scattering amplitude consists of two pieces — a 
spin-independent component proportional to XfX\ whose functional form 



a 



m a m b 



7 

(m a + m b )S - -L 



n L 7Tl r 

— i4ira — — 

<r Po 



(46) 



is identical to that of spinless scattering — together with a spin-orbit compo- 
nent proportional to 



-S b -px q 



whose functional form is 
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(47) 



We note in Eq. (|47|) the presence of an imaginary final state rescattering 
term proportional to i/po as before together with a completely new type of 
kinematic form, proportional to 1/pq which diverges at threshold. Unlike the 
term proportional to i/po, the spin-dependent piece proportional to 1/pq or 
proportional to 1/s — Sq is intriguing in that it is not imaginary but real and 
therefore does contribute to observables. (However, for bound state systems 
where the virial theorem ^mv 2 ~ - holds, this piece seems to be of the same 
order as the leading order 0(a) spin-orbit piece.) Such a kinematic form 
has been seen before by other researchers when looking at spin-dependent 
scattering. It appears in the form 



Po cos2 f Po 



\q 2 
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and has been previously identified by Feinberg and Sucher in their evaluation 
of spin-0 - spin- 1/2 scattering [TU] and by Manohar and Stewart [12]. (Note 
that in our situation, since kinematics guarantees that q 2 < Ap^ we can 
expand via 



1 1 (. q 2 q 4 



\q 2 



+ . . . 



and drop the terms higher order in q 2 jAp\ since after Fourier-transforming, 
such terms are higher order in 1/r 2 and are therefore shorter distance than 
the terms which we retain.) In any case, the presence of either of these two 
forms proportional to i/po and 1/pq would prevent us from writing down a 
well defined second order potential. 

The solution to this problem is, as before, to properly subtract the iterated 
first order potential 
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(48) 



where we now use the one-photon exchange potential ^Vq (r) given in Eq. 
(jiOj) . Splitting this lowest order potential into spin-independent and spin- 
dependent components — 
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we find that the iterated amplitude splits also into spin-independent and 
spin-dependent pieces. The leading spin-independent amplitude is 
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A^O 



and the leading spin-dependent term is 

dH (Pf 
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(52) 



which we have checked against Manohar and Stewart's expression for the 
iteration amplitude to this order [12] setting m a = m b = m. In principle we 
would also have to iterate the leading order spin-orbit piece twice. However 
this procedure yields only terms higher order in p 2 . We observe that when the 
amplitudes Eqs. (1521 and (!5T|) are subtracted from the full one loop scattering 
amplitude Eq. (jlTj) both the terms involving 1/p 2 , and those proportional to 
i/po disappear leaving behind a well-defined second order potential 
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(53) 



We observe then that the second order potential for long range Coulomb 
scattering of a spinless and a spin-1/2 particle has one component which is 
independent of the spin of particle b and which is identical to the poten- 
tial found for the case of spinless scattering, accompanied by a spin-orbit 
interaction with a new form for its classical and quantum potentials. 

Finally we would like to note that vertex corrections on the side of the 
spin-1/2 particle, particle b, give corrections to the g-factor of the spin-1/2 
particle altering the tree level value gf' = 2 to its O(a) corrected value 
= 2 + -. Since the g-factor is implicitly a parameter of the spin-orbit 
coupling piece of our leading order amplitude Eq. fl39|) and potential Eq. 
( HOI the vertex correction which we have neglected will yield a long range 
contribution of 0(a 2 ) with the same distance dependence as the leading order 
contributions proportional to L-S b /r 3 . We will neglect these contributions for 
now and include them later by writing our results for particles with arbitrary 
charges and g-factors in Appendix [Dj Therefore by using the physical values 
of the mass, charge and g-factor of the scattered particles it is sufficient to 
only consider the two-photon exchange diagrams displayed in Fig. [2] when 
considering the leading long distance corrections. 



3.2 Spin-0 - Spin-1 

It is tempting to speculate that if we extend our considerations to higher spin 
then this pattern continues — a spin-independent component identical to the 
spin-0 - spin-0 potential, accompanied by a spin-orbit interaction which is 
the same for all spins, plus additional terms which have no spin-0 or spin-1/2 
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analog. In order to test this hypothesis, we move to spin-0 - spin-1 scattering, 
and we take the spin-1 particle to be a W + boson. In order to determine the 
correct interaction vertices we must recall that the electroweak interaction is 
a non-abelian gauge theory [16]. This means that the spin one Lagrangian 
which contains the charged-W has the Proca form — 



C 



0, 



[IV) 



+ 



m 



2 



(54) 



but the SU(2) field tensor U^ v contains an additional term on account of the 
required gauge invariance [16] 



>iu, — D^Uv - DvUp - ic S u(2)Un x. U u (55) 

where csu(2) is the SU(2) electroweak coupling constant. This additional 
term in the field tensor is responsible for the interactions involving three and 
four W-bosons and for an "extra" interaction term which has the form of an 
anomalous magnetic moment and, when added to the simple Proca moment, 
increases the predicted gyromagnetic ratio from its naive value — g^± e = 1— 
to its standard model value — g^ ± = 2. As discussed in [T7] there are various 
theoretical reasons why elementary particles have a g-factor g = 2. The 
resulting one- and two-photon vertices are then found to be 
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for an incoming massive spin-1 particle with momentum p 3 and Lorentz index 
a and the outgoing one with momentum p^ and Lorentz index (3. 

We assign particle b to be the massive spin-1 particle with the incoming 
polarization vector e\ satisfying e\ ■ p 3 = and the polarization vector e b j 
satisfying ^ ■ Pa = 0. The lowest order scattering amplitude then has the 
form 
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Now we rewrite this expression using the identity 
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where we have defined the spin vector 

Sbfi 



2m, 



(59) 



The leading one-photon exchange amplitude can then be written as 
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in the symmetric center of mass frame and the transition amplitude assumes 
the form 
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The spin-independent and spin-orbit terms here are identical in form to those 
found in the spin-0 - spin- 1/2 case but now are accompanied by new terms 
which are quadrupole in nature, as can be seen from the identity 
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The corresponding lowest order potential is then 
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w : T b : s = w c T^ d s d 

and which agrees precisely with its spin-1/2 analog — Eq. fj40|) — up to tensor 
and quadrupole corrections. 

The calculation of the one loop corrections proceeds as before, but with 
increased complexity due to the unit spin. We find 
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Combining, we find the complete one loop amplitude 
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Using the identity Eq. (1581) this becomes 
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Notice here that without the e^-* ■ qe\ ■ q terms, Eq. (170]) has an identical 
structure to that of the case of spin-0 - spin-1/2 scattering — Eq. ([42]) — 
provided we substitute u(p 4 )u(p 3 ) - 



Finally, taking the nonrelativistic limit we find 
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As found in the earlier calculations, there exist terms involving both i/p 
and 1/pq which prevent the defining of a simple second order potential. The 
solution now is well known — subtraction of the iterated first order potential. 
Since the form of the spin-independent — e* B ■ — and spin-orbit — S b ■ p x q— 
terms is identical to that found for the case of spin- 1/2, it is clear that their 
subtraction goes through as before and that the corresponding pieces of the 
second order potential have the same form as found for spin-1/2. However, 
there are now two new pieces of the amplitude, the quadrupole structure 
q:T b :q which multiplies terms involving both i/p and l/p$ and the tensor 
structure p : T b : p multiplying only i/po. In order to remove these we must 
iterate the full first order potential including these quadrupole and tensor 
components. However, we find that our simple nonrelativistic iteration fails 
to remove them! We suspect the reason to be the presence of the tensor 
structure p : T b : p in the lowest order potential which is in some sense a 
relativistic correction but which when iterated yields also quadrupole pieces 
q : T b : q. A fully relativistic iteration should thus be performed which 
we will not include in this paper. It would be interesting to investigate if 
the requirement to cancel all i/po and 1/pq forms in the quadrupole and 
tensor pieces could clarify the ambiguity in the iteration of the leading order 
potential as discussed for the spinless case. 

Since we did not perform the proper iteration of the quadrupole and 
tensor pieces we merely include the spin-independent and spin-orbit pieces 
in the resulting second order potential 
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where l V^\f) denotes the tensor pieces not explicitly shown. Comparison 
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with the corresponding form of ^Vq (r) given in Eq. fl53|) confirms the 
universality which we have suggested — the spin-independent and spin-orbit 
terms have identical forms. The next task is to see whether this universality 
applies when both scattered particles carry spin. For this purpose we consider 
the case of spin- 1/2 - spin- 1/2 scattering. 



4 Spin-Dependent Scattering: Spin-Spin In- 
teraction 

In order to check universality when both scattering particles carry spin and 
to study possible spin-spin interactions, we now consider the case where both 
particles a and b carry spin-1/2. 



4.1 Spin-1/2 - Spin-1/2 

For this calculation the vertices have been given previously and the calcula- 
tion proceeds as before. The one-photon exchange amplitude is given by 
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Using the spin-1/2 identity Eq. (j34p and its pendant for particle a 
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for particle a, the one-photon exchange amplitude becomes 
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We observe that in addition to the spin-orbit pieces found previously, a spin- 
spin interaction is also present. Taking the nonrelativistic limit and working 
in the center of mass frame, we find 
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whereby the lowest order potential for spin- 1/2 - spin- 1/2 scattering is of the 
form 
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Note that since the piece proportional to S a ■ S b in Eq. ( 1761) is analytic in q 2 
it only gives a short distance contribution which is omitted in the potential 

OZZD- 

In this case when we evaluate the loop diagrams (a)-(e) of Fig. |2j we 
find that part of the spin-spin structure piece contains the form q 2 S a ■ S b 
multiplying the nonanalytic structures L and S. Due to this extra factor of 
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q 2 in this form, we must expand all loop integrals to one order higher in q 2 
than before in order to be consistent. This has been done and does make 
a difference in our results, but the very lengthy expressions for the vector 
and tensor integrals prevent us from explicitly listing them to this order in 
Appendix lAl The results for the diagrams (a)-(e) in Fig. [2] are therj§ 
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6 The results found for the spin-1/2 - spin-1/2 case are rather lengthy, and so we quote 
the results after the identities Eqs. (|34|) and (|74|) have been used. 
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where we have defined 

U a = u(p 2 )u(pi) 

and 
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with i = a,b. The sum is found then to be 
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Comparison with the result Eq. (T4"3"l) reveals again the universality which has 
been found in other cases — the forms for the scalar density and antisymmetric 
tensor components is identical and symmetric between particles a and b. 
However, there is also a new component — a spin-spin interaction. Performing 
the nonrelativistic reduction yields for the amplitude 
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Again we verify universality: the same spin-independent piece as in the pre- 
vious calculations and now two spin-orbit coupling pieces, one for the spin of 
each particle, with again the same form as found earlier. The novel spin-spin 
coupling piece consists of the last four lines of Eq. (1821) . 

Note that unlike in the case of spin-0 - spin-1 scattering where there were 
relativistic forms i b j* ■ p e\ ■ p along with the quadrupole forms e^* • q e\ ■ q in 
the leading order potential, the leading order spin- 1/2 - spin- 1/2 potential 
contains no analog relativistic term S a -pS b -p along with the spin-spin terms 
S a - qS b ■ q. We will see that now our nonrelativistic second Born iteration 
succeeds in removing all terms involving l/pjj and i/po and we find the spin- 
spin coupling piece of the 0(a 2 ) potential. 

Due to the universalities we obtained, it is clear that the iteration for the 
spin-independent piece and the spin-orbit pieces proceeds as shown before in 
the spin-0 - spin-1/2 case. As before, the second Born amplitude is 
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where we now use the one-photon exchange potential %%V<^\r) given in Eq. 
( 177|) . Splitting this lowest order potential into spin-independent, spin-orbit 
and spin-spin components — 
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we again find that the iterated amplitude splits also into spin-independent, 
spin-orbit and spin-spin pieces. As mentioned above the calculation of the 
leading spin-independent amplitude 5 1 Amp^l^g) and the leading spin-orbit 

amplitude ^Amp^_ (q) works out exactly as in the - 1/2 case, cf. Eqs. 
( l5Tj) and (|52|) . and we will not repeat it here again. 

The leading spin-spin term of the second Born iteration amplitude is new 
and we compute 
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in agreement with the corresponding terms in [T2] in the equal mass limit 
m a — m b — m - With that, the total second Born iteration amplitude becomes 
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and we observe that when this amplitude is subtracted from the full one loop 
scattering amplitude Eq. fl82|) . all terms involving and i/p disappear 
leaving behind a well-defined second order potential 
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5 Conclusions 



Above we have analyzed the electromagnetic scattering of two charged par- 
ticles having nonzero mass. In lowest order the interaction arises from one- 
photon exchange and leads at threshold to the well known Coulomb inter- 
action V(r) = a jr. Inclusion of two-photon exchange effects means adding 
the contribution from box, cross-box, triangle, and bubble diagrams, which 
have a rather complex form. The calculation can be simplified, however, 
by using ideas from effective field theory. The point is that if one is inter- 
ested only in the leading long-range behavior of the interaction, then one 
need retain only the leading nonanalytic small momentum-transfer piece of 
the scattering amplitude. Specifically, the terms which one retains are those 
which are nonanalytic and behave as either a 2 



-q 2 or a log— q . When 
Fourier transformed, the former leads to classical (/i-independent) terms in 
the potential of order a 2 /mr 2 while the latter generates quantum mechani- 
cal (ft-dependent) corrections of order a 2 h/m 2 r 3 . (Of course, there are also 
shorter range nonanalytic contributions than these that are generated by 
scattering terms of order a 2 q 2n \J —q 2 or a 2 q 2n log — q 2 . However, these pieces 
are higher order in momentum transfer and therefore lead to shorter dis- 
tance effects than those considered above and are therefore neglected in our 
discussion.) 

Specific calculations were done for particles with spin — 0, — 1/2, — 1, 
and 1/2 — 1/2 and various universalities were found. In particular, we found 
that in each case there was a spin-independent contribution of the form 
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where L = log — q 2 and S = it 2 / \/—q 2 and with S a the spin of particle a 
and Sb the spin of particle b with projections m a and rrib on the quantization 
axis. The imaginary component of the amplitude, which would not, when 
Fourier-transformed lead to a real potential, is eliminated when the iterated 
lowest order potential contribution is subtracted, leading to a well defined 
spin-independent second order potential of universal form 
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If either scattering particle carries spin then there is an additional spin-orbit 
contribution, whose form is also universal 
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where we have defined 



S a = (S a ,m af S S a ,m ai ) and S b =(S b ,m bf 
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In this case a well defined second order potential required the subtraction of 
infrared singular terms behaving as both i/p Q and 1/pq which arise from the 
iterated lowest order potential. 

In the calculation of spin-0 - spin-1 scattering we encountered new tensor 
structures including a quadrupole interaction. Unfortunately, the subtraction 
of the i/po and tensor pieces in the two-photon exchange amplitude was 
not successful with our simple nonrelativistic iteration of the leading order 
potential so that we cannot at this time give the form of the quadrupole 
component of the potential. Further work is needed to clarify this issue. The 
corrections to the spin-spin interaction have only been calculated in spin-1/2 
- spin-1/2 scattering where we found their contributions to the scattering 
amplitude and to the potential. Since we verified these forms only for a 
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single spin configuration we have not confirmed its universality which we, 
however, strongly suspect. Of course, for higher spin configurations, there 
also exist quadrupole-quadrupole interactions, spin-quadrupole interactions, 
etc. However, the calculation of such forms becomes increasingly cumbersome 
as the spin increases, and the phenomenological importance becomes smaller. 
Thus we end our calculations here. 

One point of view to interpret the universalities of the long distance com- 
ponents of the scattering amplitudes and the resulting potentials is that if 
we increase the spins of our scattered particles, all we do is to add addi- 
tional multipole moments. The spin-independent component can then be 
viewed as a monopole-monopole interaction, the spin-orbit piece as a dipole- 
monopole interaction etc. As long as we do not change the quantum numbers 
that characterize the lower multipoles (such as the charge for the monopole- 
monopole interaction or the g-factor for the monopole-dipole interaction), an 
increase in spin of the scattered particles merely adds new interactions that 
are less important at long distances. In Appendix [D] we show explicitly how 
this multipole expansion structure arises. While it is familiar from classical 
electrodynamics — i.e. at the one-photon exchange level — we are not aware 
that this has been proven for two-photon exchange processes. 

It is interesting that the same kinds of universalities of the long range 
components of the scattering amplitudes are also found in gravitational scat- 
tering [2U] and in mixed electromagnetic-gravitational scattering [21] . 

Acknowledgements 

We would like to thank John Donoghue for many clarifying discussions 
and Walter Goldberger for important comments. This work was supported in 
part by the National Science Foundation under award PHY05-53304 (BRH 
and AR) and by the the US Department of Energy under grant DE-FG-02- 
92ER40704 (AR). 

A One loop integration in EFT 

In this appendix we sketch how our results were obtained. The basic idea is 
to calculate the Feynman diagrams shown in Fig. [2J Since our calculations 
focus on long distance effects that stem from nonanalytic contributions in the 
squared momentum transfer q 2 , we only evaluate these nonanalytic pieces of 
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the one loop integrals neglecting all short distance contributions which in- 
clude the UV divergences. In practice, that means that all one loop diagrams 
where q does not run through any part of the loop can be neglected. Further- 
more, only diagrams with at least two massless propagators yield nonanalytic 
contributions, reducing the number of contributing diagrams and thus inte- 
grals further. In the end, we need four different types of integrals for our 
calculations: Bubble integrals with two massless propagators and no mas- 
sive propagator, triangle integrals with two massless propagators and one 
massive propagator and box and cross-box integrals, each with two massless 
propagators and two massive propagators. 

For simplicity we shall assume spinless scattering. Thus for diagram (a) 
of Fig. [2] — the bubble diagram — we find 

Ampm = 2! J (2^r wu? • (92) 

All vertex functions are listed in the main body of the paper, while for 
the integrals, all that is needed is their nonanalytic behavior. The exact 
expressions for the nonanalytic components of the bubble integrals read 
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Figure 3: Momentum labels for loops in triangle diagrams. On the left, the 
massive particle a with mass m a runs through the loop, whereas on the right, 
particle b with mass m& propagates in the loop. 
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where our symmetrization convention is 
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There are two distinct triangle diagrams, (b) and (c) in Fig. (2J with 
two different masses that propagate inside the loop. The momentum label 
conventions used in all triangle diagrams are shown in Fig. [3] for the two 
cases so that the expression for the amplitude for diagram (b) for example 
reads 
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In the evaluation of the integrals we use the on-shell relations 
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256m 2 V 558 m 2 
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128m 2 
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J 
I 

11V 



12 VdxuQpPa) 



J 



25 g 4 1 g 6 \ 35 
^m 4 "^^/ + 128 



1 



■m 

2 



:| 29 g 2 _ 19 g 4 1 
2\ + 64 m ^~384 m ^ + 384 m 



256m 4 V 81 m 2 324 m 4 



27g 



1 q 



q 

! 

5g 2 



I- 



35g 2 
256m : 



-J 



24m 4 V 64 m 2 



Q 



55g 4 
768m € 



128 m 4 x 
3g 4 



I- 



128m 2 
15g 4 
256m 4 



2 m 2 



128m^ 



±1 

3 m 2 

-.i 



1 + 2— + -— ) J 

m 2 6 m 4 / 



6 VOwPpPa) 

^QiiiQuPpPa) 
^(jlPvPpPa) 
PpPvPpPa 



J 



where we have defined 



1 



1 q' 
4 m 2 



(103) 



(104) 



in order to keep our notation more compact. Note that the scalar integral 
J has been expanded in the limit g 2 <C m 2 , however the expressions for the 
nonanalytic parts of the vector and tensor integrals are exact to all orders 
in g 2 when expressed in terms of the scalar integrals / and J. The triangle 
integrals listed in Eqs. (I99H103P must be used as 



{Jj J pi Jpvi Jpvpi Jf 



pupcrj 



(105) 



P=-P2, m=m a 



in diagrams where particle a (incoming momentum pi and outgoing momen- 
tum p 2 ) propagates in the loop as sketched on the left side of Fig. [31 and 
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as 



fivpcrj 



(106) 



P=+P4i ra=m\ > 



when particle b (incoming momentum p^ and outgoing momentum p 4 ) prop- 
agates though the loop with momentum labels as seen on the right hand side 
of Fig. El 

More challenging is the calculation of the box and cross-box diagrams — 
diagrams (d) and (e) in Fig. [2j For the box diagram (d) we have 



Ampp] 



d 4 k 



(2tt) 4 k 2 (k + q) 2 ((k - p 2 ) 2 - m 2 a )((k + p 4 ) 2 - m\ 



x r« (p 4 , p 4 + fc)r« (p 4 + k, p 3 ) n ^rf (p 2 , p 2 - fc)r« (p 2 - k, Pl ) . 

(107) 



The evaluation of the box integrals has been performed earlier by others with 
Ref. [22] giving a nice treatment with some of the calculational details. Un- 
fortunately, the exact expressions for the tensor integrals become extremely 
long so that we only give the form of the vector box integral. The Passarino- 
Veltman reduction of the higher tensor integrals was performed with the help 
of computer algebra, which is highly recommended. The expression for the 
scalar box integral is [22] 



K 



d 4 k 1 

(2tt) 4 k 2 ik + q) 2 {{k - p 2 ) 2 - m 2 ){(k + p A ) 2 - m\ 



16tt 2 



L 1 , 
-2-— log 

<r VA 
, 1 1 „, 

«2tt — 



A 



,s - s 



q 2 VX 



16tt 2 



q 



A — (s — s ) 
s - s ) 

s - So 



L 

-ilix — 



2m a m b 
1 



6m a m b 

s-s Q 



q 2 2^/m a m b ^s-So \ 8m a m b 



+ O((s-s ) 2 ) 

+ O((s-s ) 2 ))9(s-so) 



(108) 
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where 

A = (s - s )(4m a m b + s - s ). (109) 

Note that the expression is exact in g 2 and we only expand it in s — Sq in our 
calculations. The vector box integral is found to be 

. d*k k" 



+ 



(2tt) 4 k 2 (k + q) 2 ((k-p 2 ) 2 -ml)((k+p 4 ) 2 - m|) 
(2m 2 + 2m a m b + s-s )J a + (2m 2 + 2m a m b + s-s )J b -AK M 
2 [((m« + m 6 ) 2 + s - s )g 2 + A] ? 
(4m a m fc +2(5-5 )+g 2 )J a + (4mg-g 2 )J fe +g 2 (2mg+2m a m b +g-g )A" ^ 

2[((m a + m 6 ) 2 + s-s )g 2 + A] 
(4m 2 -g 2 )J a +(4m a m i ,+2(s-s )+g 2 )J6+g 2 (2m 2 +2m a m b +s-s )A' 



P2 



2[((m a + m b ) 2 + s - s )q 2 + A] 



-Pa 
110) 



with Jj = 

g 2 , s - s 



J\ , and we notice that its denominator vanishes in the limit 

\m=m,i 7 

— > 0. More specifically, the denominator can be written as 



2 [{{m a + m b ) 2 + s - s )q 2 + A] 
2 



m l + m l + 2 



8Po 



8pg 



m 2 + m 2 + 2 



m 2 + m 2 + 2 



Pl + 



pI + 



m 2 +P 2 , 



m 2 +p 2 



M - q 2 ) 



1 — sin 2 



cos 



(111) 



where we have used Eq. ( fTOj) and g = — q = — 4pgsin |. We see that 
the denominator vanishes for p —>■ and for backward scattering at = ir. 
Unless we consider backward scattering where the denominator vanishes and 
thus the amplitude diverges, we have 4p^ > q 2 , and since p\ originates from 
the relativistic structure s — sq, we therefore expand our vector and tensor 
box integrals first in q 2 and then in s — Sq. Denominators that vanish in 
the limit g 2 , s — Sq — > are a common feature for all box vector and tensor 



-s ) 



integrals, and they are the source of the l/(s 
results for the scattering amplitude. 

In the case of the cross-box diagram (e) the amplitude reads 

f d 4 k 1 

~ J (2^ k 2 (k + q) 2 {{k - p 2 ) 2 - m 2 )((k - p 3 ) 2 - m 2 ) 



I/Pq components in our 
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x r« (p 4 , p 3 - k)rW (pg - fc, P3 ) T^yV^ (pa, P2 - fc)r« (pa -A;,pi). 

(112) 

Now we need the cross-box scalar integral which can be deduced from the 
result for the box scalar integral by replacing the set of Mandelstam variables 
(s, t) by (u, t) where t = q 2 and s + t + u = 2m 2 , + 2m 2 . Again, we only give 
the scalar and vector integrals because the exact expressions for the higher 
tensor integrals become very long. The resulting expressions are 



K'-- 



(2vr) 4 k 2 (k + q) 2 ((k - p 2 ) 2 - ml)((k - p 3 ) 2 - m 2 ) 



16tt 2 



L 



2m a m b 




s 



Q 



6m a mb 
2 ■' 2(s 



+ 



6m a m b 
?>Qm 2 a m 2 



d 4 k 



lAm a rrib 
+ 0{q 6 ) (l + O(s-s ) 



S,i +0((s- s ) 2 ) 

So) +0((s- s ) 2 ) 
+ O((s-s ) 2 ) 



5m a m b 
9(s-s ) 



+ 



(2tt) 4 k 2 (k + q) 2 ((k - p 2 ) 2 - m 2 )((k - p 3 ) 2 - m 2 ) 
{2m 2 a -2m a m b -(s-So)-q 2 )J a + (2ml-2m a m b -(s-So)-q 2 )J b -AK' ^ 
2[((m a + m b ) 2 + s-So)q 2 + A] Q 
(4m a m b +2(s-So)+q 2 )J a -(4m 2 h -q 2 )J b -q 2 (2m 2 b -2m a m b -(s-So)-q 2 )K' 

2 [((m a + m b ) 2 + s- s )q 2 + A] 
(4m 2 a -q 2 )J a -(4m a m b +2(s-So)+q 2 )J b + q 2 (2m 2 a -2m a m b -(s-s )-q 2 )K' 



P-2 



2[((m a + m b ) 2 + s-So)q 2 + A} 



-Pi 

;ii3) 



where 

A = (s - s + q 2 )(4m a m b + s - s + q 2 ). (114) 

We point out that we did not include an imaginary component in the case 
of the cross-box scalar integral whereas for the box integral in Eq. (11081) 
we included both an imaginary and a real part. The reason for that is that 
the ^-function multiplying the imaginary part in Eq. (llOSp for the cross-box 
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integral becomes 6{u — s ) and it vanishes in the kinematic region we are 
considering since 

u = {m a - m b f - (s - s ) - q 2 < s = (m a + m h f . 

In this way all amplitudes quoted in the text can be generated. 

B Fourier Transformations 

In this appendix we collect all Fourier transformation integrals needed to 
evaluate the potentials in coordinate space. 

d 3 q 1 1 



(2tt) 3 |gf Anr 
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(2tt) 3 \q\ ' 7T 2 V r 6 r 4 

A -iq-r, M 2 1 

e q log \q\ 
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(2tt) 3 ^ 27rr5 
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(27r) 3 ' ! ' 7rr 
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C Iteration Integrals 

In this appendix we evaluate the integrals 
[H; H r ; H r ,] = i j ^- e * illlM^L « (116 ) 

which are needed in order to perform the iteration of the lowest order Cou- 
lomb potentials. Note that we have introduced a small photon mass A 2 -C p$ 
in order to avoid singularities, but since we are only interested in the long 
distance effects we do not show the singularities in A in our expressions. The 
evaluation of the integral H has been given by Dalitz as [23J 

tf = ,W^^ = -W^4. (117) 

Po q po q 

In order to determine the vector integral H r we define 

H r = A{ Pi +p f ) r (118) 
and contracting with (p, +Pf) r , we find 

2A(fi + p f ) 2 = (2A 2 + Ap 2 )H - Am r Y - X( Pi ) - X(p f ) (119) 

where 

Y = - / 77T^ = ~T = -2a 2 5 120 



and 

,2 



X (p f ) = X( P/ ) = le 



dH e 2 



(2ir)-|f_fl|. + A »A-£ + ie 

m r iX , 

= -247ra 2 ^log 121 
Po 2p + «A 

Note that the integrals X only depend on pi or pj and therefore do not yield 
any terms nonanalytic in q 2 . Thus we drop the contributions of the X's and 
we have 



A = -—^-^[AplH-Am r Y]^a 2 (^S- l 2^^\ (122 
8pg(l-sm 2 f) L J \p 2 Po q 2 ) 
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In the case of the tensor integral we define 

H rs = B5 rs + C (pi +Pf)r(j>i +Pf)s + D (pi ~Pf)r(Pi ~Pf)s (123) 

and we require three conditions in order to evaluate the coefficients B, C and 
D. Neglecting again the integrals X, these are 

i) 

S rs H rs : 3B + (Apl - f )C + q 2 D ~ p 2 H - 2m r Y 



ii) 



(pi + p f ) r H rs : B + {Apl - q 2 )C ~ — ^- [p 2 // - m r Y] - m r Y 

1 4^ 



iii) 

Solving, we find 



{ Pi -p f ) r H rs : B + q 2 D~0 



B * -£(h-?£ 
4 V Po 



D - ^tf-^H (124) 



4 V" Vl 

Keeping only the leading terms in q 2 we have then 

tt ■ a 2 m r L 

H ~ -lAna - 

iJ r ~ (Pi + p f ) r -«27ra + a —^S + . . . 

_ 2 / . 2 m r L 1 2 m r 

-Hrs - q l^OL ^ - -OL —^b + . . . 

V Po Q 2 



+ (Pi +P/)r(Pi +P/)s -«7TQ! + —=-5 + . 

V Po 1 Po 

+ {Pi~ Pf)r(Pi~ Pf) s ( -ina 2 — ^ + \a 2 —^S + . . .) (125) 

V p q 2 2 p 2 Q J 
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D Generalized Results and Interpretation 



How can we interpret the universalities we have found? As an example, let 
us first consider two spinless charged particles of charge e where the leading 
order Coulomb interaction between these two charges is 

V(r) = -. (126) 
r 

Now if we replace one of the two charges by a spin- 1/2 particle of charge — 2e 
the spin-independent leading order Coulomb interaction becomes 

—In 

V{f) = —. (127) 
r 

We see that the universality of even the leading order Coulomb potential de- 
pends on having equal charges. In this section, we extend our calculations to 
arbitrary charges and g-factors, and our results lead us to the interpretation 
that the universalities originate from a multipole expansion of the long range 
scattering amplitudes and potentials. 

D.l One-photon Exchange Potential 

It has long been known that a particle with spin 5* has 2S + 1 multipole 
moments j2l] . The one-photon exchange potential thus exhibits a multipole 
expansion as we know it from classical electrodynamics: 

Spin-0 — Spin-0 

The Lagrangian for a spin-0 particle with arbitrary charge q = Ze reads 

C = (iD^iDy - m 2 f (128) 
with = <9 M + ieZA^ and the Feynman rules for the vertices become 

^ 1) (P2,Pi) = -iZe(p 2 +pi)f, 

T$fa,Pi) = 2z(Ze)\ l/ . (129) 

The one-photon exchange potential for a spin-0 Particle a with mass m a , 
charge q a = Z a e and a spin-0 particle b with mass m b , charge q b = Z b e is 
then 

V«(f) ~ (130) 
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and exhibits merely a monopole-monopole interaction, Coulomb's law, pro- 
portional to 1/r. 



Spin-0 — Spin-1/2 

Now we introduce the Lagrangian for a spin-1/2 particle of arbitrary charge 
q — Ze and arbitrary g-f actor g, 



L = i,{iip- m)i> - - F^a^ (131) 



with again = <9 M + ieZA^, which yields the Feynman vertex rules 

i m . Ze(g — 2) 

2 ^y{P2,Pi) = -iZe-i lx + — o w {v2-V\Y 

H^Vi) = 0. (132) 



The one-photon exchange potential for a spin-0 particle a with mass m a , 
charge q a = Z a e, g- factor g a and a spin-1/2 particle b with mass m b , charge 
q b = Z h e, g-factor g b is 

Vc [T) ~ —T-XfXi-—^ ^1 L-S b (133) 

where L- S b = (r x p) ■ S b = r ■ (p x S b ). Thus, besides the leading monopole- 
monopole interaction - Coulomb's law - proportional to 1/r, we observe an 
additional monopole-dipole interaction, the spin-orbit coupling. The coeffi- 
cient of the monopole-monopole interaction does not depend on the g-factor 
whereas the spin-orbit piece does. 

Using the physical values for the masses, charges and g-factors then in- 
cludes all one-photon exchange effects to all orders in a, i.e. if we take 
g = 2 + ^ (for a particle of charge ±e), we take into account the 0(a 2 ) 
long-distance contribution from the one loop vertex correction diagram. 



Spin-0 - Spin-1 

For the spin-1 case we consider a particle of charge Ze with g-factor g in the 
Lagrangian 

£ = -\ul u U^ - m 2 ^ + tZe(g - l)0j>„F^ (134) 
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where = — D v (p^ with = + ieZA^. Once the charge and 

the g-factor are determined, the quadrupole moment from this Lagrangian 
is fixed. If one wants to include an arbitrary quadrupole moment one has to 
add a dimension 6 operator [25] which we will not do here since it complicates 
the following two-photon exchange calculations considerably. The resulting 
Feynman rules read 

lr Sa(P2,Pl) = iZe Va/3(P2+Pl)n 

-v^(gp2 -(g- i)pi) a - Vfj,a(gpi -(g- 1)22)0 

vJ/9a(P2,Pi) = -i{Ze) 2 [2r]^r] a p - rj^p - rj^T] va \ (135) 

and the one-photon exchange potential for a spin-0 particle a with mass m a , 
charge q a = Z a e and a spin-1 particle b with mass 777.6, charge q b = Zbe, 
g-factor g b is 

i T/ (i) M Z a Z h a ^ ^ Z a Z b a (g b - l)m a + g b m b -> - 

+ ^3(2, b -l) f:T&:F 
r 5 4m^ 

where we neglected relativistic terms involving e^* • p q • p. Now besides 
the monopole-monopole and monopole-dipole pieces seen before, a new piece 
of monopole-quadrupole structure constitutes the highest multipole in the 
expansion for the spin-1 particle. 

Spin-1/2 - Spin-1/2 

The one-photon exchange potential for a spin-1/2 particle a with mass m a , 
charge q a = Z a e, g-factor g a and a spin-1/2 particle b with mass m b , charge 
q b = Z b e, g-factor g b is 

rQ-)/x\ ~ Z a Z b a 



Z a Z b a (g b - l)m a + g b m b r - _ a] a 



r 

.n ( riu — 1 lrr?,_ 4- n^m.i, 

,, -L-S b XfXi 



«a"«b 



Z a Z b a g a m a + (g a - l)m b ^ - 6f b 
r 3 2m\m b a ^ ^* 
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Z a Z b a g a g b 
r 5 4m a m b 



(3S a ■ fS b ■ f - r 2 S a ■ S b 



.) (137) 



In this case, we observe a monopole-monopole piece, two monopole-dipole 
pieces aka spin-orbit pieces and a dipole-dipole piece, the spin-spin interac- 
tion. 

Clearly, the one-photon exchange potentials (and scattering amplitudes) 
for two charged particles of various spins exhibit a multipole expansion, and 
as one would expect for a multipole expansion, higher spins only add higher 
multipole interactions while all lower multipole interactions are universal, i.e. 
of identical form as for lower spins. That then implies that the numerical 
coefficients in the multipole expansion do not depend on structures as for 
example S 2 or on coefficients that characterize higher multipoles, for example 
Coulomb's law cannot depend on the g-factors but only on the charges and 
the spin-orbit interaction does not depend on the quadrupole moment. 

D.2 Two-photon Exchange Potential 

At the two-photon exchange level the amplitudes and potentials we calculated 
exhibit the same universalities as found in the one-photon exchange case 
where they are explained in terms of a multipole expansion. For particles 
with arbitrary charges and g-factors the results for the second order potentials 
read 
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+ 



(Z a Z b a) 2 ({g b -2)ml + (2g b - 3)m 2 a m b + 2(g b - l)m a m 2 + g b m^j 
2m 2 a ml(m a + m b )r A 

(Z a Z ba ) 2 h((-3g 2 + lQg b -18)m a + (-3g 2 + 16g b -A)m b 



8itm'lm b i r b 



L-S b 



+ (f) 



Uv^ 2 \r) 
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a] a b] b 

"X.f Xi "X.f "X.% 
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Am 2 a m\[m a + m b )r 4 
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(Z a Z b a) 2 S a ■ fS h ■ f/r 2 



+ (g 2 a + 20g a -12)g b m 2 b 
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8m 2 a m 2 h {m a + m b )r 4 
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+ 9 2 a 9b + ^9a9b(9a + 9b) ~ 56g a g b + Kol + 9b) + %a + 9b) + 16 

5(Z a Z b a) 2 hS a ■ f S b ■ r /r 2 
x — - — ^— — (138) 

where we have introduced the short notations SO a = L- S a XfX\ an d SOb = 
xf%fL-S b . 

The two-photon exchange potential exhibits a similar structure as the one- 
photon exchange potential, a "generalized multipole expansion" . Its classical 
part of O(h ) starts off with a monopole-monopole piece proportional to 
1/r 2 , then there is a monopole-dipole piece proportional to L ■ S/r A etc. 
Thus the "generalized multipole expansion" of the classical part is similar 
to the multipole expansion of the one-photon exchange potential, but it has 
one additional power of r in the denominator. The "generalized multipole 
expansion" of the quantum O(h) part of the two-photon exchange potential 
however is seen to start with a monopole-monopole piece that falls off as 1/ r 3 
followed by monopole-dipole pieces that go as L ■ S/r 5 etc. 

That then suggests the interpretation of the universalities we have found 
for the long-distance two-photon scattering potentials and amplitudes as fol- 
lowing from a "generalized multipole expansion" where "generalized" means 
that the multipole expansion of the potential does not start with a monopole- 
monopole term proportional to 1/r as for a usual multipole expansion but 
proportional to l/r n with n > 1. 

It would be interesting to see if one could prove this multipole expan- 
sion scheme and thus the universalities found using low-energy theorems for 
Compton scattering amplitudes and combining two Compton scattering am- 
plitudes to a two-photon exchange scattering amplitude using dispersion rela- 
tions. Moreover, one could speculate that a three-photon exchange potential 
would exhibit a similar structure with a "generalized multipole expansion" 
and universalities. 
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